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ABSTRACT 

The  free  lateral  vibration  of  a nonlinear  viscoelastic  beam-column 
subjected  to  an  initial  compressive  axial  load  is  considered. 

The  constitutive  law  is  formulated  with  a linear  elastic  term  and 
with  power  functions  of  stress  in  the  transient  and  steady  creep 
terms,  and  is  of  the  nonlinear  generalized  Kelvin  type. 

. 

By  assuming  that  the  stress  caused  by  the  oscillation  is  of  much 
smaller  magnitude  than  the  initial  stress,  the  problem  is  linearized. 
The  problem  is  analyzed  for  five  special  viscoelastic  models  using 
small  deformation  theory,  and  numerical  results  are  discussed  for 
a stainless  steel  alloy. 


1.  INTRODUCTION 


As  noted  by  Hoff  creep  in  metals  must  have  a damping  effect  on 
mechanical  vibration,  since  it  absorbs  energy.  If  one  is  considering 
a vibration  problem  in  a temperature  range  for  which  creep  may  be 
present,  the  damping  due  to  creep  should  be  included  in  the  analysis. 
Cozzarelli,  Wu  and  Tang(^)  considered  a beam  constructed  of  a ma- 
terial which  may  be  represented  by  a nonlinear  Maxwell-Kelvin  model, 
and  which  was  subjected  to  a free  lateral  vibration  while  under  initial 
axial  tension.  In  this  paper  we  extend  their  analysis  in  two  ways; 
the  first  of  which  is  to  employ  a more  general  viscoelastic  law  of  the 
nonlinear  generalized  Kelvin  type  in  order  to  provide  greater  flexibi- 
lity in  the  fitting  of  experimental  data.  And  secondly  we  extend  the  type 
of  loading  to  include  the  case  of  initial  axial  compression.  The  growth 
or  decay  of  damped  vibration  under  compression  is  a problem  of  prac- 
tical importance  in  studies  of  structural  reliability  and  safety,  and 
thus  we  give  particular  attention  to  this  aspect  of  our  study.  Creep 
collapse  of  a fuel  element  in  a nuclear  reactor  is  an  example  of  such 
problems  which  is  of  current  interest. 

The  particular  case  considered  here  is  r uoar  viscoelastic  beam- 

column  which  is  simply  supported  and  unu^.  a initial  axial  compres- 
sive load.  The  beam-column  is  assumed  to  be  initially  straight  and  of 
an  isothermal  homogeneous  medium  with  a rectangular  cross-section. 
By  assuming  that  the  initial  stress  increment  is  of  much  higher  mag- 
nitude than  the  stress  increment  caused  by  the  oscillation  of  the  beam, 
the  problem  is  linearized  in  a manner  analogous  to  that  in  ref.  2. 

The  incremental  approach  to  the  theory  of  dynamic  creep  stability  in 
columns  has  been  developed  in  a number  of  papers;  of  these  we  men- 
tion in  particular  Rabotnov  and  Shesterikov(^»  ^),  Jahsman  and  Field^, 
Jahsman'  ',  and  Dist^fano  and  Sackman^’  The  present  development 
more  closely  parallels  that  in  ref.  7-8,  since  it  is  based  on  a nonlinear 
viscoelastic  relation  rather  than  on  a mechanical  equation  of  state  as  in 
ref.  3-6.  Our  specific  goal  is  to  demonstrate  in  detail  how  for  actual 
values  of  the  creep  parameters  one  may  readily  construct  the  complete 
incremental  displacement  solution  for  rather  general  constitutive  rela- 
tions primarily  by  means  of  simple  hand  calculations. 

The  governing  displacement  equation  is  derived  in  Section  2 using 
small  deformation  theory,  and  the  solution  of  this  equation  for  a non- 
linear Maxwell-Kelvin-Kelvin  (or  simpler)  material  is  discussed  in 
Section  3.  This  letter  section  employs  the  property  of  physically  small 
damping  due  to  creep  to  develop  simple  but  accurate  estimates  to  the 
roots  of  the  characteristic  equations.  Numerical  results  for  a stainless 
steel  are  considered  in  Section  4,  where  we  examine  the  effects  of 
steady  creep,  slow  transient  creep  and  rapid  transient  creep. 
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2.  DERIVATION  OF  THE  GOVERNING  EQUATION 


2.  1 Constitutive  Daw 


A one -dimensional  constitutive  law  of  the  nonlinear  generalized  Kelvin 
type,  similar  to  the  law  discussed  in  detail  and  used  in  ref.  9,  will  be 
employed  here  in  its  integral  form 
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Here,  T,  £(T)  and  C(T)  represent  time,  strain  and  stress  respectively 
and  E is  the  modulus  of  elasticity;  X is  the  steady  creep  parameter  and 
n is  the  steady  creep  power;  Z ^ are  retardation  times,  are  transient 
creep  parameters  and  are  transient  creep  powers;  and  sgnC^T)  is  the 
usual  signum  function.  By  setting  N=1  in  eq.  (1)  we  arrive  at  the  consti- 
tutive equation  used  in  ref.  2,  except  that  here  the  use  of  the  signum 
function  enables  us  to  consider  negative  stresses. 

Consider  a material  with  an  initial  constant  stress  °rQ  (positive  or  nega- 
tive), subjected  to  a small  disturbance  at  T=TQ.  This  disturbance  is 
translated  into  a small  additional  increment  of  stress,  i.  e. 


cr(T)  = <ro  + £^(T)  (2) 

where  (T ) is  the  stress  increment  and  £ is  a small  quantity. 

It  then  follows  that 
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Also,  as  a result  of  the  disturbance  the  strain  assumes  the  form 

£(T)=  £o(T)+  £6  (T)  + . . . . (4) 

Substituting  eqs.  (2)-(4)  into  eq.  (2),  subtracting  the  relation  for  the 
initial  state,  and  retaining  terms  of  order  L , we  obtain  the  following 
linear  constitutive  relation  in  the  increments: 
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Beam -Column 

As  a specific  example  we  consider  the  vibration  of  an  initially  straight 
Euler-Bernoulli  beam-column,  which  is  subjected  at  T=0  to  a pres- 
cribed axial  force  Pq  (tensile  or  compressive)  followed  later  by  a 
small  lateral  disturbance  at  T=TQ.  Using  the  perturbation  technique 
employed  in  obtaining  eq.  (6),  the  classical  strain  displacement  and 
equilibrium  relations  for  such  a beam -column  yield  the  following  rela- 
tions in  the  first  order  terms: 
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In  eqs.  (7),  X is  the  neutral  axis  of  the  beam-column  which  vibrates 
in  the  XY  plane,  V (X,  T)  is  the  lateral  displacement,  M = -fjr  YdA  is 
the  bending  moment  for  a beam-column  with  cross -section  A,  p is  the 
linear  density,  and  W is  the  lateral  load  which  will  be  used  to  deter- 
mine the  initial  conditions  of  the  free  vibration. 

For  convenience,  we  introduce  the  nondimens ional  quantities 

v T-T 

v X o T ... 

v = — ± . x = r » t=  r~  » w = — r w (a) 

V*  L T*  pv* 

Here,  V*  is  some  convenient  constant  reference  displacement,  L is 
the  length  of  the  beam,  T is  a convenient  factor  which  has  the  dimen- 
sion of  time,  and  v,  x,  t and  w are  the  non-dimensional  lateral  dis- 
placement, axial  coordinate,  time  and  lateral  load,  respectively. 

Combining  eq  i.  (6),  (7)  and  (8)  we  finally  obtain  the  governing  equation 
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and  where  the  second  moment  of  area  is  I = S^{  dA.  Note  that  the  non- 
dimensional  force  a has  the  same  sign  as  PQ , whereas  b,  c,  d^  and  j- 
are  all  positive  nondimensional  material  parameters.  Also  note  that  as 
a result  of  material  nonlinearity  the  parameters  c and  d^  depend  on  I P I . 

In  Section  3,  eq.  (9)  will  be  analyzed  for  five  different  models.  These 
models  are,  in  order  of  decreasing  complexity,  nonlinear  Maxwell- 
Kelvin-Kelvin  C *°  (dj  and  — * 0)  for  i >,  3)  , nonlinear  Maxwell- 
Kelvin  and  jj  ■*  0)  for  i > 2 ] , nonlinear  standard  solid 

£ (d^  and  j.-»  0)  for  i V 2 and  X-*«*  (c  -*■  0)]  , nonlinear  Maxwell 

[ f. -»  M (d^  and  j.  -*  0)  for  i * 1 ) and  linear  elastic  L E’.  -*  «•  (d^  and  -*  0) 
for  i > 1 and  >>  -*«>(c  -*  0)]  . The  differential  form  of  eq.  ^9)  for  the  non- 
linear Maxwell-Kelvin-Kelvin  (M-K-K)  model  is  obtained  via  successive 
differentiation  as 
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The  remaining  differential  forms  for  the  other  special  models  are 
easily  obtained  from  the  limits  listed  above. 
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Next,  the  nondimcnsional  loading  function  is  expressed  in  the  form 

w - s(*)  n (t)  ( 


where  g(x)  is  a function  of  x , ami  •?  (t)  is  a doublet  function  which 
excites  a free  vibration  for  t ^ 0 . For  the  solution  of  eq.  (11)  with 
forcing  function  (12),  it  is  convenient  to  set  w=0  in  eq.  (11'  and  use 
initial  conditions  at  t=0+.  For  these  conditions  we  integrate  eq.  (11) 
successively  between  t=0-  (where  all  quantities  vanish)  and  t=0^  and 
obtain 
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(!3  c) 
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We  see  that  the  beam-column  has  in  effect  been  deformed  laterally 
t=0+  in  the  form  of  g(x)  and  then  released  from  rest. 

3.  SOLUTION  OF  THE  GOVERNING  EQUATION 

3.  1 The  Characteristic  Equation 

As  previously  noted,  we  shall  consider  the  free  vibration  of  a simply- 
supported  nonlinear  Maxwell-Kelvin-Kelvin  (M-K-K)  beam-column 
(and  simpler).  Noting  the  boundary  conditions 


v(0,t)  = v(l,t)  = 


v(0,  t) 


v(l,  t) 


and  the  fact  that  all  the  space  derivatives  in  governing  eq.  (11)  are  of 
even  order,  we  shall  take  g(x)  of  eq.  (12)  in  the  convenient  form 

g(x)  = sin  (mix)  (1  5) 

where  m indicates  the  mode  of  vibration.  For  m = l we  see  from  eq.  (13  a) 
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and  (15)  that  v(l/2,0+)  = 1,  which  in  turn  requires  lor  this  case  that 
V*  in  eq.  (8)  equal  the  initial  lateral  displacement  at  mid-span. 


A solution  to  eq.  (11)  then  follows  simply  as 
v(x , t)  = T(t)  sin  (m T x) 


(16) 


where  T(t)  is  a function  of  time.  Substituting  eq.  (16)  into  eq.  (11_) 
with  the  load  w=0,  we  obtain  an  ordinary  differential  equation  in  T(t)  as 
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The  characteristic  equation  for  the  nonlinear  M-K-K  material  is  ob- 
tained by  substituting 


T (t)  = 


(19) 


into  eq.  (17).  From  this  point  on  we  shall  confine  our  attention  to  the 
case  of  initial  compression  (a  = - |a|  ),  for  which  the  characteristic 
equation  is  given  by 

b^5  + b(c  + d:  + d2  + jL  + j2)$4 
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where  the  o(.  are  the  six  coefficients  in  this  fifth  order  equation.  The 
characteristic  equation  for  the  simpler  models  are  obtained  by  allow- 
ing the  appropriate  creep  parameters  to  tend  toward  zero,  as  described 
in  Section  2.2.  For  convenience,  the  coefficients  in  the  characteris- 


tic  equations  for  all  of  the  models  are  summarized  in  Table  1.  Note 
that  in  the  summation 


p-i 


the  quantity  p is  the  order  of  the  equation,  which  equals  2,  3,  3,  4 and 
5 for  the  E,  M,  S-S,  M-K  and  M-K-K  materials  respectively. 

Having  obtained  the  roots  to  the  various  characteristic  equations,  the 
general  solution  to  eq.  (11)  for  unrepeated  roots  then  follows  from  eqs. 
(16)  and  (1  9)  as 

v(x,  t)  = / i c.e  sin(m7Tx)  (21) 

i=l  1 


the  c^  are  constants  which  may  be  determined  from  initial  con- 
-ions,  eq.  (13). 
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3.2  Discussion  of  Roots  ,-tnd  Solutions  for  Initial  Compression 

In  this  section  we  examine'thc  nature  of  the  roots  to  the  five  charac- 
teristic equations  listed  in  Table  1,  and  also  discuss  the  correspond- 
ing solutions  in  accordance  with  eq.  (21).  beginning  with  the  clastic 
material  we  see  that  the  two  roots  are  imaginary  for  |al<i  and  real 
(one  positive)  for  |a|>  1.  Thus  the  solution  for  la!  > 1 exhibits  expo- 
nential growth,  cind  we  recognize  the  value  lal=l  [i.  c.  |PQl  = m^T^El/l^ 
v*a  eq.  (10)  and  eq.  (18)}  as  the  critical  buckling  load  for  an  clastic 
column.  On  the  other  hand,  the  solution  for  |a|<l  oscillates  with  the 
nondime  ns  ional  natural  frequency 


which  decreases  as  the  initial  load  increases.  Since  actual  structures 
are  usually  designed  to  operate  at  loads  below  this  critical  load,  we 
shall  assume  that  for  all  the  models 


III  < 1 (23) 

Now  we  consider  as  a group  the  four  other  viscoelastic  models  listed 
in  Table  1,  and  remark  that  the  damping  due  to  creep  is  physically 
small  (see  ref.  2).  Accordingly,  the  solutions  for  all  these  viscoelastic 
models  will  contain  two  terms  which  oscillate  at  essentially  the  elastic 
frequency  to  as  given  by  eq.  (22)  and  decay  slowly  with  time.  Such  terms 
correspond  to  two  roots  with  negative  real  parts 
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1,2 


-p  + ico 


where  to  »p> 


(24) 


which  implies  that  the  creep  parameters  are  of  small  order  /3  when 
compared  with  to.  If  we  substitute  eq.  (24)  with  (22)  into  the  character- 
istic equation  for  the  nonlinear  M-K-K  material  and  ignore  terms  of 
higher  order,  we  find  that  the  imaginary  part  is  identically  satisfied 
while  the  real  part  yields  the  relation 


/» 


C * dl  * d2 

2(1  - III) 


(25) 


Eq.  (25)  applies  equally  well  to  the  nonlinear  M-K  material  (d.,=0), 
the  nonlinear  S-S  material  (c=d_=0),  and  to  the  nonlinear  M material 

(dl=d2=°). 

The  remaining  p-2  roots  will  now  be  examined  for  the  various  models 
with  the  aid  of  Descartes  " rule  of  signs.  We  see  from  Table  1 that  in 
all  four  models  the  coefficients  ocq  and  are  positive,  and  is  also 
positive  due  to  inequality  (23).  For  the  Maxwell  material  we  note  the 
**3  is  always  negative,  and  thus  the  remaining  third  root  for  this  mate- 


L- 
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rial  must  be  positive  real,  and  the  total  solution  contains  two  decay- 
ing osciLatory  terms  and  one  exponentially  growing  term.  Wc  may 
obtain  a good  estimate  for  this  positive  real  root  by  solving  the  cha- 
racteristic equation  with  lj  = 0,  since  this  is  equivalent  to  dropping 
the  inertia  and  the  associated  oscillatory  terms.  In  so  doing  we  obtain 


t * laJjE. 

^3  1 -la  l 

For  the  standard  solid  material  we  see  that  (X^  may  he  positive  or 
negative,  depending  on  the  value  of  I a \ . Thus,  if 


(26) 


\al  > a 


where  ac  = J1/(d1+j1)<  1 


(27) 


the  remaining  third  root  must  be  positive  real,  whereas  if  | a|  < a 
this  root  is  negative  real.  The  estimate  to  this  r-  ot,  as  obtained  by 
setting  1j=0,  is  given  by 


|at  (dx  + j1)-j1 


1 - I al 


(28) 


which  exhibits  the  above  mentioned  behaviour  in  regard  to  sign.  In  con- 
trast to  the  M material,  the  exponential  term  in  the  total  solution  for  the 
S-S  material  grows  in  time  only  if  the  load  is  large  enough  for  inequality 
(27)  to  be  satisfied.  Thus  for  this  material  there  is  a second  critical 
load  |P  | defined  by  )a|  = ac,  which  upon  using  equations  (10),  (18)  and 
(27)  yields  the  nonlinear  expression  in  |P  l 


lPocl  L 

m/l 


2 _ 


*1  1 Poc/A< 

V1 


qi-i 
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+ E 


= 1 


(29) 


The  quantity  in  brackets  is  the  effective  reciprocal  elastic  modulus  of 
the  linear  spring  in  series  with  the  nonlinear  Kelvin  element  (with  2^=0). 
Note  that  lPocl  is  less  than  the  critical  elastic  buckling  load,  and  also 
that  due  to  material  nonlinearity  the  elastic  modulus  of  the  Kelvin  ele- 
ment depends  on  |Poc)  itself. 

Consider  next  the  nonlinear  Maxwell  -Kelvin  material,  where  there  are 
two  roots  in  addition  to  the  complex  pair  given  by  eq.  (24).  Again  exa- 
mining Table  1,  we  see  that  whereas  cx^  is  always  negative,  may  be 
either  positive  or  negative.  However,  in  both  cases  there  is  only  one 
variation  in  sign,  and  thus  there  is  one  positive  real  and  one  negative 
real  root.  Setting  b=0  in  Table  1,  we  obtain  the  following  estimate  to 
these  two  real  roots  \ ft 

-j1+|a|(c+d1+j1)2+  [[  jj  -|a|(c+d1+j1)]2+  4(1  -ISlJlal^c} 

^3,4  “ 2(1  - | al  ) 


$4  < 0 


(30) 
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which  simplifies  to  cq.  (26)  and  cq.  (2  8)  when  the  appropriate  creep 
parameters  are  sen  equal  to  zero. 


Finally,  we  consider  the  nonlinear  M.ixwell-Kelvin -Kelvin  material, 
for  which  there  are  three  roots  in  addition  to  the  two  previously  dis- 
cussed complex  roots.  If  0^  in  Table  1 is  negative,  then  for  the  cor- 
responding values  of  |a|  one  can  prove  that  (X^  must  also  be  negative. 
However,  if  is  positive,  then  <X^  may  be  cither  positive  or  negative. 
For  these  three  cases  there  is  still  only  one  sign  variation  in  the  coef- 
ficient sequence,  and  thus  there  is  again  one  and  only  one  positive  real 
root.  The  other  two  roots  may  in  principle  be  either  both  negative  real 
or  a second  complex  conjugate  pair,  but  our  numerical  calculations 
have  yielded  no  characteristic  frequency  other  than  the  one  near  the 
elastic  frequency  (see  Section  4).  The  three  real  roots  are  estimated 
by  the  solutions  to  the  cubic  equation 


3 

K + 


r,  . <c+di+d^l2 

f (j1c+j2c+d1j2+d2j1)|ar| 

|_ji+j2 - i - isi  _p 

+ Ljlj2'  1 - lal  J 

Jlj2c|a| 


1 - lal 


= 0 


(31) 


U> o 


An  examination  of  the  test  function  for  Cardan's  formula  confirmed  the 
observation  that,  for  typical  values  of  the  parameters,  the  three  solu- 
tions to  eq.  (31)  are  in  fact  real.  Thus,  the  overall  behaviours  of  the 
M-K-K,  M-K,  S-S  (with  |al  >ac)  and  M materials  are  essentially  the 
same,  i.  e.  there  is  an  initial  period  of  decaying  oscillation  which  is  fol- 
lowed by  a period  of  exponential  growth. 


For  convenience,  the  character  and  the  estimates  of  the  roots  for  the 
various  models  are  summarized  in  Table  2.  The  actual  evaluation  of 
the  coefficients  c^  in  the  solution  [eq.  (21)^  from  these  roots  is  some- 
what tedious  for  the  more  general  models,  and  thus  we  will  obtain  the 
results  for  the  M and  S-S  models  only.  Utilizing  Table  2 we  write  the 
displacement  solution  for  both  these  models  as 


r .At  £3*1 

v(x  , t)  = sin  mVx  (c^cos  u)t  + c 2 sin  cot)  e + c^e 


(32) 


Also,  initial  conditions  (13  a-c)  with  (15)  and  (22)  yield 

2v 

at 


a2 

v(x,  0 ) = sin  mTpc , ~ (x,  0 ) = 0,  ^-(x,  0+)  = sin  m7rx  (33) 


dt 


Evaluating  c^,  c^,  c^  in  eq.  (32)  from  conditions  (33),  we  obtain  the 
complete  solution  for  the  M and  S-S  models  as 
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/(x.,  t)=5 


sin  mTfx 
2.  2 


(f  3+(*r+« 


( ^ 3 + 2/i  53+  w ) cosut  + 


A(?2+  ^A+fe2) 


co 


sinwt 


. c + ft  c 


,-A‘  . 


(34) 


Note  that  although  the  coefficient  C3  is  very  small  (since  /%«. co)  it  is 
always  positive,  and  thus  for  positive  £3  the  displacement  for  t large 
is  in  the  same  direction  as  the  initial  displacement. 


For  linear  viscoelastic  materials  it  is  well  known  that  the  S-S  material 
is  characterized  by  two  critical  loads,  and  that  for  materials  with  a 
long-time  viscous  behaviour  (e.  g.  M,  M-K,  M-K-K)  the  lateral  displa- 
cement always  grows  with  time  (e.  g.  see  Distefano,  ref.  10).  Accord- 
ingly,  the  results  obtained  here  concerning  the  roots  for  a particular 
nonlinear  material  are  consistent  with  previously  established  results 
from  the  theory  of  linear  viscoelasticity.  The  roots  which  we  have  es- 
timated in  this  section  will  be  obtained  by  numerical  computer  calcu- 
lations in  the  next  section  for  a stainless  steel  beam  at  two  levels  of 
load  and  temperature. 


4.  NUMERICAL  RESULTS 


In  order  to  illustrate  the  roots  and  solutions  discussed  in  the  previous 
section,  we  consider  a 36  in.  long  stainless  steel  (type  316)  beam  with 
a rectangular  cross-section  which  is  2 in.  wide  by  3 in.  deep.  For  the 
axial  force  PQ  and  the  temperature  0 we  consider  two  cases,  i.  e. 

PQ  = -39  000  lb  with  9=  1 300°F  and  PQ  = -60,000  lb  with  6 = 1 500°F, 
where  we  shall  soon  see  that  the  latter  case  satisfies  | a|  > ac  Qsee 
eq.  (27)3  while  the  former  satisfies  I a | < ac.  The  material  properties 
p,  n,  qj  , as  listed  in  Table  3 (with  lb,  hr  and  in.  units)  are 

identical  with  those  obtained  in  ref.  2 for  a M-K  model  using  the  data 
of  Garofalo  et  al.  ^).  The  values  of  E listed  in  Table  3 for  the  two 
temperatures  were  obtained  by  extrapolation  from  data  at  lower  tem- 
peratures given  in  ref.  13,  and  they  are  slightly  less  than  those  used 
in  ref.  2. 


It  remains  for  us  to  obtain  some  values  for  q£,  yU  ^ , i.  e.  the  creep 
parameters  in  the  second  Kelvin  element  of  the  M-K-K  model.  Reexa- 
mining the  data  in  ref.  11,  12  we  see  that  the  relationship  between  the 
"initial"  strain  and  the  stress  is  clearly  nonlinear  and  greater  than  the 
calculated  linear  elastic  strain.  We  assume  here  that  the  difference 
between  this  initial  strain  and  the  calculated  linear  elastic  strain  is  a 
time -dependent  strain  due  to  rapid  transient  creep.  In  the  absence  of 
more  detailed  data  we  set  q£  = q^  and  ^ x 10"^,  and  then  calcu- 

late from  the  data  in  ref.  11,12  at  stresses  less  than  10,000  psi. 
Although  these  estimated  values  arc  rather  crude,  they  serve  the  pur- 
pose of  enabling  us  to  illustrate  the  effect  of  rapid  transient  creep.  All 
material  parameters  for  both  temperatures  are  listed  in  Table  3. 


TABLE  3 - Material  Properties 


0 (°F) 

1300 

1500 

p,  10'10  lb-hr2/in2 

3.45 

1 

3.45 

E,  107  lb/in2 

2.0 

1.5 

n = q]  = q2 

3.64 

3.  50 

X,  104  lb/in2 -hr1/11 

19.18 

7.  86 

fx,  , 104  lb/in2 

3.  10 

2.  11 

04 

a 

•H 

\ 

jQ 

O 

< 

7.45 

5.26 

o 

X 

II 

r—i 

333 

21.7 

TABLE  4 - Nondimensional  Constants 


j2=jl  x 10 
a_ 


P = -30,  000  lb 
©°  = 1 300°F 

P = -60.  000  lb 
9°  = 1 500°F 

4.  377  x 10-2 

1.167  x 10"1 

8.  566  x 10  ^ 

1. 142  x 10'6 

6.  943  x 10“6 

1.071  x 10‘3 

1.  585  x 10“5 

4.  925  x 10'3 

6.  516  x 10“3 

2.  013 

8. 342  x 10'3 

1.280  x 10"1 

5.006  x 10'2 

2. 592  x 10‘3 
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I or  the  evaluation  of  the  nondimensional  material  constants  from 
eqs.  (10)  and  (18)  ve  select  m = l (first  mode)  and  T*  = 1 /3600  hr, 
which  converts  real  time_  scale  T-T0  from  hours  to  seconds  feq.  (8)]. 
Material  constants  lal,  b,  c,  dj,  d£,  jj,  plus  the  critical  load  a 
^eq.  (2 7 ) "3  arc  listed  in  Table  4 for  the  two  previously  mentioned  cases 
of  load  and  temperature.  Note  that  for  PQ  = -30,  000  lb  at  6=  1300°F  we 
have  lal  < ac<  1,  whereas  for  PQ  =-60,  000  lb  at  0 = 1 500°F  we  have 
t a I ^ 1 • Thus,  while  both  cases  are  below  the  critical  linear  elas- 
tic load  as  required  by  eq.  (23),  only  the  former  case  is  also  below 
the  additional  critical  load  for  a standard  solid  material. 

Using  the  nondimensional  constants  listed  in  Table  4,  the  roots  to  the 
characteristic  equations  given  in  Table  1 were  obtained  by  a program 
on  a high  speed  computer.  The  roots  for  PQ  = -30,  000  lb  at  6 = 1300°F 
are  presented  in  Table  5,  and  the  roots  for  PQ  = -60,  000  lb  at  0 = 1 500°F 
are  given  in  Table  6.  When  these  numerically  computed  roots  are  com- 
pared with  the  simple  estimates  summarized  in  Table  2,  we  find  that  the 
estimates  are  remarkably  accurate  and  in  fact  agree  with  the  computer 
obtained  roots  to  four  significant  figures.  Finally,  as  an  illustrative  case 
we  give  below  the  complete  numerical  solution  for  the  M model  at 
Pc,  = -60,000  lb  with  6 = 1 500°F  as  obtained  from  eq.  (34)  and  Table  6: 

v(x.t)  = sinTx{4.  755  x 1(T13  exp(1.416  x 10’4  t)  + exp(-6.  064  x 104t) 

• [cos(879.  4 t)  + 6.  896  x 10~7  sin  (879.  4 t)l  V . (35) 

The  deflection  solutions  based  on  the  roots  listed  in  Tables  5 and  6 are 
in  complete  conformity  with  the  predictions  of  Section  3.2.  As  expected, 
all  solutions  contain  two  terms  which  oscillate  at  the  elastic  frequency  cu 
["eq.  (22)]  , although  the  sin  cot  term  is  negligible  in  the  present  exam- 
ple since  the  motion  starts  from  rest  C e.  g.  see  eq.  (3  5)]  . As  predic- 
ted by  eq.  (25)  for  0 , the  rate  of  decay  of  these  oscillatory  terms  in- 
creases as  the  load  and  temperature  increases  and  as  additional  creep 
components  are  added  to  the  model.  The  material  with  the  larger  va- 
lues of  | PQ|  and©  (Table  6)  and  with  the  rapid  transient  creep  terms  in- 
cluded (i.  c„  M-K-K)  exhibits  sufficient  damping  (0=1.  143)  to  cause 
complete  decay  of  the  oscillatory  terms  over  a period  of  several  seconds. 

Examining  the  real  roots  in  Tables  a end  6,  we  see  that  as  predicted 
only  the  S-S  material  at  the  lower  values  of  |PQ|  and  0 (Table  5)  does 
not  have  an  associated  positive  real  root.  In  this  one  case  the  perturbed 
displacement  tends  toward  zero  as  and  thus  this  material  under 

these  conditions  may  be  termed  "stable".  In  the  other  seven  cases  the 
perturbed  displacement  tends  toward  infinity  as  t-*oo,  and  accordingly 
these  may  be  considered  "unstable".  One  can  show  that  the  coefficient 
of  the  exponentially  growing  term  is  always  positive  as  in  eq.  (3  5),  and 
thus  in  these  cases  the  column  "collapses"  (after  a period  of  decaying 
oscillation)  in  the  same  direction  as  the  initial  displacement. 
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TABLE  5 - Roots  for  the  Various  Models  at 
PQ  = -30,  000  lb  6 = 1 300°F 


Mate  - 

rial 

?li2  = -f!  ±i“ 

?3 

*4 

^5 

ft 

to 

E 

- 

1.057xl03 

- 

- 

- 

M 

3.  630x1 0'^ 

1 . 057x1 03 

3. 178xl0"7 

- 

- 

S-S 

8.  283xl0“6 

1 . 057xl03 

-1.  085xl0-7 

- 

- 

M-K 

1.  192xl0"5 

1 . 057xl03 

6.301xl0"7 

- 

M-K-K 

3.450xl0'3 

1.057x103 

6. 628x1 0'^ 

-8.  117xl0~3 

-4. 148xl0"7 

TABLE  6 - Roots  for  the  Various  Models  at 

P =-60,000  lb  0=  1 500°F 
o 


Mate  - 
rial 


E 

- 

8. 794x1 02 

- 

M 

6.064x10  4 

8.  794xl02 

1. 416xl0~4 

S-S 

2.  788x10  3 

8.  794xl02 

6. 380xl0-4 

M-K 

3.  394x1 0 3 

8.  794x1 02 

7. 81 9x1 0"4 

M-K-K 

1.  143 

8.  794x1 02 

1. 396x1 0'1 

5.  CONCLUDING  REMARKS 


We  have  demonstrated  in  detail  how  one  may  use  simple  hand  calculations  to 
construct  the  complete  incremental  displacement  solution  for  a vibrating 
nonlinear  viscoelastic  beam  under  initial  axial  compression.  Solutions 
were  obtained  for  rather  general  constitutive  relations,  including  those 
with  the  long-time  viscous  behaviour  which  is  characteristic  of  metals 
at  high  temperatures.  For  such  cases  the  lateral  displacement  always 
grows  with  time,  even  if  the  load  is  less  than  the  critical  load.  The 
solutions  obtained  may  be  used  to  estimate  a critical  time  at  which  the 
lateral  displacement  exceeds  a prescribed  acceptable  value. 
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APPENDIX  - LIST  OF  SYMBOLS 


a,  b,  c,  d.,  j 
a,  b 
a 

c 


A 

c. 

1 

E 

£ 

g{  x) 

I 

L 

m 

M(X,T) 

n 

N 

P 

P 

o 

P 

oc 


t 

T 


W(X,  T) 
w(x,  t) 
x 
X 

Y,  Z 


nondimensional  material  constants,  eqs.  (10) 

nondimens ional  matcrml  constants  for  mode  m,  eqs.  (18) 

nondimensional  critical  load  for  standard  solid  mate- 
rial, eq.  (27) 

cross-sectional  area  of  beam 
constants  of  integration,  eq.  (21) 
modulus  of  elasticity 
small  quantity 

space  dependence  of  nondimensional  loading  function, 
eq.  (12) 

second  moment  of  area 

length  of  beam 

vibration  mode  number 

bending  moment 

creep  power,  eq.  (1) 

number  of  Kelvin  elements 

order  of  characteristic  equation.  Table  1 

initial  axial  force  on  beam 

critical  load  for  standard  solid  material,  eq.  (29) 
creep  power  of  i-th  nonlinear  Kelvin  element,  eq.  (1) 
nondimensional  time,  eqs.  (8) 
time 

constant  reference  time,  eqs.  (8) 
time  at  beginning  of  disturbance 
function  of  time 

nondimensional  lateral  displacement,  eqs.  (8) 
lateral  displacement 

constant  reference  displacement,  eqs.  (8) 
lateral  load  on  the  beam 
nondimensional  lateral  loading,  eqs.  (8) 
nondimensional  axial  coordinate,  eqs.  (8) 
centroidal  axis  of  beam 
coordinates  in  cross-section  of  beam 
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List  of  Symbols  (continued) 


magnitude  of  real  part  of  complex  roots,  cqs.  (24) 
and  (2  5) 

* 

€ 

strain 

lit) 

doublet  function 

» 

6 

tempe  rature 

X 

steady  creep  parameter,  eq.  (1) 

Mi 

transient  creep  parameter  of  i-th  nonlinear  Kelvin 
element,  eq.  (1) 

K 

root  of  characteristic  equations,  Table  1 

e 

linear  density 

Cr 

stress 

* 

r. 

l 

retardation  time  of  i-th  Kelvin  eiement,  eq.  (1) 

4 

CO 

nondime  ns  ional  frequency,  eq.  (22) 

( )o 

indicates  initial  value  during  interval  T < T 

o 

• 

C) 

indicates  increment  due  to  lateral  disturbance 

n 

indicates  d/dt 
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